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Abstract 



A study of Final State Effects (FSE) on the dynamic structure function 
of superfluid ^He in the Gersch-Rodriguez formalism is presented. The main 
ingredients needed in the calculation are the momentum distribution and the 
semidiagonal two-body density matrix. The influence of these ground state 
quantities on the FSE is analyzed. A variational form of p2 is used, even 
though simpler forms turn out to give accurate results if properly chosen. 
Comparison to the experimental response at high momentum transfer is per- 
formed. The predicted response is quite sensitive to slight variations on the 
value of the condensate fraction, the best agreement with experiment being 
obtained with uq = 0.082. Sum rules of the FSE broadening function are also 
derived and commented. Finally, it is shown that Gersch-Rodriguez theory 
produces results as accurate as those coming from other more recent FSE 
theories. 
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I. INTRODUCTION 



Deep inelastic neutron scattering (DINS) has been extensively applied to the study of 
quantum fluids, since Hohenberg and Platzman'Jll proposal of using DINS to determine 
the momentum distribution n{k) of helium atoms in superfluid ^He. The determination 
of n{k) in quantum liquids is a challenging problem of fundamental interest .1 In fact, the 
knowledge of n{k) provides very useful information to understand basic properties of the 
quantum nature of these systems as the Bose-Einstein condensation. At the same time, 
the theoretical analysis of DINS probes and stimulates the development of modern many- 
body techniques. These issues have been the main motivations of a considerable amount of 
measurements and theoretical work on liquid ^He and other quantum liquids.il"0 

The inelastic scattering of neutrons by liquid ^He is described by the double differential 
cross section 



where b is the scattering length, and kf are the initial and flnal wave vectors of the 
scattered neutron , and q and uj are the momentum and energy transferred from the neutron 
to the sample. The dynamics of the sample is entirely contained in S{q,u!), the dynamic 
structure factor, which is the Fourier transform of the density-density correlation function.^ 
At sufficiently high momentum transfer, the scattering is entirely due to single atoms and 
S{q,u!) can be accurately described by the Impulse Approximation (IA),0 provided that the 
interatomic potential does not have an infinite repulsive core. In this regime, the kinetic 
energy of an atom recoiling from a neutron collision is much larger than the potential energy 
due to the interaction with the neighboring atoms, so that collisions of the former with other 
atoms can be neglected. The lA predicts a simple relation between S{q,u!) and n{k), 



where ur = q^ /2m is the free atom recoil frequency, m is the mass of the '^He atoms and n{k) 
is the thermally averaged occupation probability of the single particle state of momentum 
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S{q,uj) , 



(1) 





k, which reduces to that of the ground state at T = 0. The delta function in Eq.(0) takes 
care of the momentum and energy conservation in the scattering event between the neutron 
and a single atom. Assuming S{q,uj) = Sia{(1,uj), the momentum distribution n{k) can be 
extracted from Eq. (|^) by simple differentiation. Notice that, in the previous equation and 
henceforth, h is set to 1. 

In isotropic systems, where n{k) depends only on the modulus of k, it is useful to 
introduce the Compton profile 

JiA{Y) = -^SjA{q,u;). (3) 
m 

which is driven by a single variable 

Y = ^(u-^], (4) 

and fulfills Y scaling.El If a finite fraction of atoms uq occupies the zero momentum state, 
Jia{Y) presents a 6 peak of strength no at F = 0. However, this expected signature of 
the condensate, is not observed in experiments performed at momentum transfer as high as 
23A~^,i because the lA spectrum is broadened by both final state effects (FSE) and instru- 
mental resolution effects (IRE). Hence, the theoretical interpretation of the experimental 
data requires not only the knowledge of n{k), but also an accurate description of both the 
dynamics which determines FSE and the instrumental broadening function.! 

Several methods to account for FSE have been proposed.lii'@Hll Among them, we will 
focus on the so-called convolutive theories, in which 

/oo 
duj'SiA{q,uj')R{q,uj - uj'), (5) 
-oo 

where R{q, uj) is the FSE broadening function. 

After the first attempt^ to approximate R{q, uj) by a Lorentzian of width proportional 
to the ^He-^He cross section, Gersch et ahi expressed the response function S{q,iu) in a 
1/q series expansion, whose coefficients are given by integrals of many-body correlation 
functions averaged on the ground state of the system. In this approach, the response when 



g — > oo is given by the first term of tlie expansion of the incoherent part of S{q,uj), which 
turns to be exactly the lA. However, the theory could not deal with realistic interatomic 
potentials presenting a strong repulsion at short distances. To overcome this problem, 
Gersch and Rodrigue:^! proposed a cumulant expansion of S{q, t) which provides an adequate 
frame for calculating the response function at high momentum transfer. The full calculation 
is impractical, but with some approximations based on physical grounds, S{q, u) can be 
expressed in terms of the one- and the semidiagonal two-body density matrices, and the 
two body interaction. At the time the theory was proposed the numerical application was 
made with a very simple approximation of the two-body density matrix that resulted in an 
overestimation of the response at the quasielastic pealcS 

The main purpose of the present work is to revisit Gersch-Rodriguez theory, and show 
that using a realistic two-body density matrix one gets a S{q, u) in good agreement with 
both experimental data and more recent theories of FSE.Ill0 

In the next section, a review of the theory is presented. Section III is devoted to the 
discussion of the results and their comparison with the experimental data. A sum rules 
analysis of R{q, u) is presented in section IV. In Section V our results are compared with 
other FSE theories, and finally section VI summarizes the main conclusions of the work. 



II. GERSCH RODRIGUEZ THEORY OF FSE 

In the Gersch-Rodriguez theory,! the density-density correlation factor S{q, t) is ex- 
pressed as the product of the lA and the FSE correcting function by means of a cumulant 
expansion. The n-th order cumulant accounts for the correlations among the struck atom 
and clusters of n particles in the background. In the high momentum transfer limit, those 
terms with n = 1 carry the most significant corrections. At this level, the FSE broaden- 
ing function can be expressed as a function of the interatomic potential and the one- and 
two-body density matrices. 

The starting point in Gersch-Rodriguez theory is the time representation of the response 
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NS{q,t) = Y.< e-'^^^'e'^^e'^^^^e-'"^ >=Y.< e^'»('"^-'"')e-''i'"^e^^*e''i'"^e-^^* > , (6) 

3,1 3,1 



which can be brought to the following form 



NS{q, t) = e'"^"' < e^''^'"^-'"'^ e'^^' T exp i f dt' H{rj - v/)) 



(7) 



where T is the time-ordering operator and H{rj — Vqt') is the actual Hamiltonian of the 
system where the position coordinate of particle j has been shifted by an amount v^t'. As 
the interatomic potential considered is velocity-independent, one can write 



if(r,-v/) = i/ + f/,(v/, , 



(8) 



with 



and 



3 «<J 



(9) 



Uj,mivgt') = [V{rj - Vgt', Vm) - V{rj, r^)] 



(10) 



where Vq = q/m and ujg = q^ /2m. 

The incoherent part of the response, which is defined by taking particles labeled j and 
/ in Eq.(0) to be the same, is the only contribution at large q. In this limit, S{q,t) may be 
written in the following way 



S{q,t) = e'"'* < e'"^*P^e'^* Texp 



-'^ m^l 



(11) 



where tJ{\qt') is the previously defined potential operator but with the position operators 
evaluated at time t' rather than at t = 0. Notice that expression ([TlD is as hard to evaluate 
as the original S{q, t). An exact treatment would require the knowledge of the time evolution 
of the whole system, so different approximations should be made in order to deal with this 
last relation. 



Gersch and Rodriguez! performed a cumulant expansion of the ground state expectation 
value of Eq.(|TT]). The expansion contains an infinite number of terms, and allows for the 
factorization of the lA from the total response 



S{q,t) = SiA{q,t) R{q,t) 

R{q, t) being the FSE correcting function given by 

1 



(12) 



R{q, t) = exp 



E< 



JtVqPl 



1 -Texp z / dt'Ui^m{^,t'] 



> +... 



(13) 



Up to this point, the result is exact because it is nothing more than a rearrangement 
of the different terms entering in S{q,t). The first problem in the calculation of Eq.(|13|) is 
associated to the infinite number of terms appearing in the exponential. Such a difficulty can 
be skipped if one looks for the underlying physics contained in each term: the contribution of 
the n-th order cumulant to S{q, t) accounts for the correlations between n-particle clusters 
during their interactions with the struck atom. One may expect that the first significant 
correction to the lA is produced by the multiple scattering of the struck particle with the 
atoms of the media, considering them independently of each other. This corresponds to a 
truncation of the series beyond the first order cumulant. 

The second problem lies on the evaluation of the time-dependence appearing in the 
particle coordinates of f/i,m(V(,t'). In the large q limit, the displacement of the struck particle 
is much larger than the average movement of the background atoms. Thus, one can discard 
the time dependence of r(t) in Ui^m- This is a safe procedure as, even though the inclusion 
of such a time dependence avoids hard-core collisions between the struck particle and other 
target atoms, the contribution to R{q, t) coming from those situations vanishes due to rapid 
oscillations in the imaginary exponential of Eq.(|r^). Therefore, one can writei 

—j— f drp2{r, 0; r + v^t, 0) [l - exp {t f dt' {V{v + v,(t - t')) - V{v + v,t))l 

Ol[Vqt) J L I JO J , 

(14) 



R{q, t) = exp 



where pi and p2 are the one-body and semidiagonal two-body density matrices!^ of the 



system, respectively. R{q, t) is a complex function, but its Fourier transform is real because 
its real part is even and its imaginary part odd under the change t — > —t. 

Equation (0) predicts S{q, t) as the product of Sia{(1, t) and R{q, t), and therefore S{q, u) 
is the convolution of SiA{q,^^) and R{q,uj) 

/oo 
duo' SiA{q-,^^')R{qi^ — i^') ■ (15) 
-oo 

In the particular case of liquid ^He, the momentum distribution n{k) may be written as 

n{k) = (27r) Vo <5(k) + n(A;) , (16) 

where no is the condensate fraction value and n{k) stands for the occupation of non-zero 
momentum states. Consequently, Sja{(1,^) is split in two parts 

Cq^ \ 777, roo / qr2 \ _ 

^^^2^)^ A^q %-|| ^ ""^^^ '^^^'''^\;'~2^j^ • ^^^^ 

where the first term on the rhs is the condensate response which appears as a delta peak of 
strength no located at the quasielastic energy, and S'/yi(g, cu) is the non-condensate contri- 
bution of n{k) to the lA. Introducing the West variable Y = mu/q — q/2, Sja{(1,^^) can be 
expressed in terms of the Compton Profile 

±SjA{q,u;) = JjAiY)=no6iY) + -^ H kn{k) dk, (18) 
m Att^p J\y\ 

which scales in Y 

Moreover, at high q the response is usually expressed in terms of Y through the relation 

J{q,Y) = ^S{q,uj) , (19) 
and thus, Eq.(^) is transformed into 

/oo roc _ 

dY'JjA{Y')R{q, Y-Y')= noi?(g, Y) + / dY'JjAiY') R{q, Y-Y') , (20) 
-oo J — oo 



where 



R{q,Y) = ^R{q,u) . (21) 
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III. NUMERICAL RESULTS 



In this section, we present results for the FSE correcting function R{q^ Y) and the re- 
sponse function J(g, Y) calculated in the framework of the Gersch-Rodriguez formalism. 
The input density matrices pi(r) and p2(i"i, r2; r'l, r2) used to calculate Jia{Y) and R{q^Y) 
have been obtained in the framework of the HNC theory0'0ii from a variational many body 
wave function containing two- and three-body correlations.ii The variational minimization 
has been performed for the HFDHE2 Aziz potential at the experimental equilibrium den- 
sity (p = 0.365cr^^, cr = 2.556A). The ground state description obtained with this wave 
function is in good agreement with recent Green's function Monte Carlo calculations.ii@ 
The discussion is separated in two parts, the first one being devoted to the study of both 
Y) and J(g, Y) and their comparison to experimental data, and the second one to the 
analysis of the dependence of these functions on the different approximations used in the 
variational description of the ground state wave function. 

The actual calculation of the FSE broadening function is initially performed in time 
representation (|T3|). R{q,x) is a complex quantity which can be written in the following 
way: 

R{q, x) = e'^(«'^) [cos{ip{q, x)) + i sin(^(g, x))] (22) 

with 



4){q,x) 



2 r 1 /"^ 

-— / (irp2(r, 0; r + x) sin^ - — / du {V (r + - u) - V (r + x.)} 

Pi[xj J [2Vq Jo 

If \ ^ 

ip(q, x) = — --^ / drp2(r, 0; r + x) sin — du {V(r -|- x — u) — V(r + x)} 

Pi (x) J v„ Jo 



(23) 



X being Vgt. As can be seen form Eq.(^), (t){q,x) and i/j{q,x) are even and odd functions 
of X, respectively. Therefore, the real and imaginary parts of R{q, x) are respectively even 
and odd under the change x —x, and consequently R{q, Y) is real. Even if the potential 
becomes very repulsive at short distances, as is the case of the Aziz potential, Eq. ( [T^ ) gives 
an R{q, Y) which does not diverge. 
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The real and imaginary parts of R{q,x) are shown in Fig. 1 for q = 23.lA~^. In 
the relevant range of x, Rei?(g, x) has a dominant decreasing behavior. The Rei?(g, x) 
and lmR{q,x) are related to the symmetric and antisymmetric components of R{q,Y), 
respectively. As the imaginary part is much smaller than the real part, R{q, Y) is mostly 
symmetric around Y = 0. In Fig. 2, we show ip{q,x) and (f){q,x) at g = 23. lA^^. (f)(q,x) 
is a negative and a monotonously decreasing function of x, causing both the real and the 
imaginary parts of R{q,x) tend to zero when a; ^ oo (p^). 

In Fig. 3, we show R{q,Y) at two different values of g, 23.1 A^^ and 15.0 A^^. The 
main trends of R{q, Y) in all FSE convolution theories are the same: a dominant central 
peak and small oscillating tails which vanish as |y| increases. As one can see, the shape of 
R{q, Y) smoothly changes with q, this variation being reflected in an overall redistribution 
of the strength between the main peak and the wings. When q increases, the peak appears 
higher and narrower pointing to the tendency of R{q, Y) to become a delta distribution in 
the limit q oo. 

The existence of a finite condensate fraction uq in superfluid ^He plays an important role 
in the FSE corrections, as is reflected in Fig. 4 where the broadening of the condensate and 
non-condensate parts of Jja{Y) are separately shown. The small differences between JiaO^) 
(dotted line) and the convolution of Jia{Y) with R{q, Y) (long-dashed line) reveal small FSE 
on the non-condensate part of the response at high q. In contrast, the broadening of the 
condensate term (short-dashed line), i.e., the convolution product of R{q,Y) and nQ6(Y), 
contributes to J{q, Y) as hq R{q, Y) which is a function with an appreciable width and 
height. The inclusion of the latter term produces a total J{q, Y) (solid line) which manifests 
a sizeable departure from the lA prediction. Therefore, FSE corrections in superfluid ^He 
appear to be relevant even at so high g's.lii 

A direct comparison between theoretical and experimental dynamic structure factors is 
not possible due to the presence of instrumental resolution effects (IRE) in the experimental 
data acquisition process. It would be desirable, from a theoretical viewpoint, to remove 
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the IRE inherent to the measured response, especially at high q where they become larger. 
However, the latter is an ill-posed problem due to the statistical noise of the data, and thus 
the only way to compare theory and experiment is by convoluting the theoretical J(g, Y) 
with an instrumental resolution function J(g, F). At present, I{q,Y) is obtained from a 
Monte Carlo simulation of the experimental setup, and in contrast to earlier models used 
in neutron scattering analysis, it is neither Gaussian nor symmetric around F = 0, and is 
comparable in width and height to -R(g, Y) at those momenta.0 The influence of /(g, Y) in 
the response is sketched in Fig. 5 for q = 23.1 A^^. As one can see, the introduction of 
the IRE in the response (solid line) appreciably modifies J{q,Y) (dashed line). The most 
important effect of I{q, Y) is to quench the central peak reducing the effects of the FSE 
correction on Jja(Y), whereas the tails remain almost unchanged. 

In Fig. 6, we present results of J(g, Y) broadened by the IRE at different values of 
q in comparison with inelastic scattering data at T = 0.34 K from Ref. ^. There is an 
overall agreement between the predicted and the observed scattering data, the quality of 
the Gersch-Rodriguez theory being comparable to results provided by other theorieJil0 
(see also section V). It is worth to notice that all FSE theories are stressed when applied to 
intermediate q values. This is also apparent in our results, as one can see for the lowest q 
value reported in Fig. 6. Thus, whereas the experimental peak shifts its location to a small 
negative Y value, the theoretical one is shifted to so small positive values of Y that it is not 
appreciable in the figure. 

The most relevant quantity in the calculation of J{q, Y) is the momentum distribution 
n{k) which completely determines the Compton profile Jia{Y). The influence of n{k) in 
J(g, Y) is shown in Fig. 7. The dashed and solid lines correspond to a Jastrow n{k) {nj{k)) 
and a Jastrow plus triplet n{k) {njT{k)), respectively. The condensate fraction predicted by 
the two approximations are slightly different, Uq = 0.091 and = 0.082. This reduction 
of uq, produces a small decrease of strength in the peak of J{q, Y) bringing our theoretical 
prediction closer to the experiment. A basic ingredient in the calculation of R{q, Y) is the 
semidiagonal two-body density matrix, which in the framework of the HNC theory is given 
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P2(ri, r2; r'l, = ppi{rn')gwd{ri2)9wd{ri'2) exp [A{ri, r2; r'l)] 



(24) 



where piirw) is the one-body density matrix, gwd{f) is an auxihary two-body radial distri- 
bution function and A{ri, r2; r'l) is the sum of the Abe diagrams. Notice that the structure 
of p2 allows for the exact cancellation of pi in Eq. (|14|). As the explicit dependence of p2 in 
no is introduced in pi, the influence of no in R{q,Y) is almost negligible. We have verified 
that the inclusion of three-body correlations does not appreciable modify the structure of 
R{q, Y). Consequently, three-body correlations can be omitted in the calculation of R{q, Y). 
In a further step, we have also studied the influence of the Abe diagrams using a Jastrow 
wave function. As is well known, it is not possible to calculate y4(ri, r2, r'^) exactly but a 
good estimation of its contribution can be obtained through the scaling approximation .El 
The inclusion of the Abe diagrams in Eq. (p4D using the scaling approximation produce neg- 
ligible effects in the final form of R{q, Y). In fact, the Abe terms, which quickly vanish when 
the inter-particle separation increases, only modify the structure of p2 when coordinates 1, 
1' and 2 are very close to each other. These small changes in p2 are then suppressed when 
integrated to obtain R{q,t). Furthermore, one can slightly change the functions gwd{f) and 
no influences in R{q, Y) are observed. This fact, which will be explicitly commented in Sec. 
V, points to the relevance of the functional decomposition of p2 rather than the exact form 
of the functions entering in it. 

IV. SUM RULES 

In this section we study the sum rules satisfied by the Gersch-Rodriguez FSE broadening 
function R{q,Y). From the relation 



and the first sum rules of both Jia{Y) and the incoherent part of J{q,Y), an equivalent 
set of y-weighted integrals for R{q, Y) can be derived.^ Notice that equation ( p5|) can be 




(25) 
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taken as a possible definition of R{q., Y) provided tliat q is large enough for the coherent 
part of J(g, Y) to be negligible. These sum rules are model independent, and so any suitable 
convolutive FSE broadening function must fulfill them. The first sum rules of -R(g, Y) are 

/CXD 
dYR{q,Y) = l 

/oo 
dYYR{q,Y) = 
-oo 

/oo 
dYY^R{q,Y) = 
-CXD 

/CO iryi P 

^ dYY'R{q, Y) = —pj dvg{r){ci ■ V) V(r) . (26) 

As we are only considering the incoherent part of the response, mQ{q) is 1 at any q. Both 
the first and second moments of R{q, Y) vanish because the Impulse Approximation exactly 
fulfills the incoherent sum rules. Finally, the third moment of R{q, Y) is expressed in terms 
of the two-body radial distribution function g{r) and the interatomic potential V{r), which 
are not included in Jja{Y). 

Relations (|26|) are exact and partially define the behavior of R{q, Y). Therefore, one can 
use them to check the accuracy of R{q, Y) calculated using different approximations. In 
the Gersch-Rodriguez theory, the sum rules analysis can be analytically performed. In fact, 
expressions for the sum rules can be easily derived from the time derivatives of R{q, t) at 
t = 



1 d 



dt^ 



(27) 



Performing a McLaurin expansion of R{q,t), Eq. (0), the different coefficients of the 
series are directly related to the F-weighted sum rules. In this way, one obtains the relations 

R,GR/ \ 1 

mi (g) = 
m^'^^(g) = 

mf^'^'^iq) = ^ / rfrp2(r, 0; r)(q ■ V) V(r) + ^ / rfr [(q ■ V)F(r)] [(q ■ V.)p2(r, 0; x)]^^ (2^ 

(J P J (J P J 
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where rrij^ ' (g) stand for the F-weighted integral of the FSE function in Gersch-Rodriguez 
theory. Integrating by parts the second term of m^'^^[q), and taking into account general 
symmetry properties of p2, one can express ra^''^^{q) in the following way 

mf''''{q) = ^/ rfrp2(r,0;r)(q- V)V(r) . (29) 

As the diagonal part of p2 is p'^g{r), the analytic expression of m^{q) is recovered. Therefore, 
the zero, first, second and third moments of R{q, Y) are exactly fulfilled in the Gersch- 
Rodriguez theory. 

Nevertheless, the exact p2 is not known, and the use of an approximation can produce 
numerical differences between Eqs. (pOf ) and (^61). In fact, we have checked that the inclusion 
of the Abe terms in the variational p2 defined in Eq. (|^) is crucial in reproducing g{r) in 
its diagonal part, and consequently mf{q). 

V. COMPARISON WITH OTHER FSE THEORIES 

FSE theories can be classified in different groups depending on the way they incorporate 
the corrections to the lA. The two most important groups are, on one hand, convolutive 
theories in which the total response is expressed as a convolution of J/a(^) and R{q,Y) 
and, on the other, additive theories where the leading FSE corrections are summed up to 
the lA. Examples of theories belonging to the first class are those of SilveJi or Carraro and 
Koonin.0 An example of additive theory is that originally derived by Gersch, Rodriguez 
and Smith,! which was next generalized by Rinatlli to treat also hard core potentials. 

Gersch-Rodriguez formalism was the first in predicting convolutive FSE corrections. 
Silver's and Carraro and Koonin's theories appeared some years after. In this section, we 
present a comparison between their results and our predictions obtained in the framework 
of the Gersch-Rodriguez theory. 

In the Gersch-Rodriguez theory R{q, Y) is formulated in terms of the semi-diagonal 
two-body density matrix of the system. In the present work, a variational ansatz for this 
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quantity has been employed and discussed, but at the time the formahsm was developed 
only a qualitative description of p2 "was available. This led the original authors to use a form 
of p2 based on a Hartree-Fock approximation and the Schwartz inequalityi 



P2(ri,r2;r'i,r2) = ppi{riv)^ g{ri2)g{rv2) , (30) 

Pi(r) being the one-body density matrix and g{r) the two-body radial distribution function. 
At that time, detailed microscopic calculations of g{r) were not available, so they had to 
approximate it. The form selected for the radial distribution function was simply a step 
function 

g{r) = e{r - ro) , (31) 

with a parameter tq to mimic the radius of the hole of g{r). Originally, tq was taken as a fit- 
ting parameter. However, theoretical arguments brought them to fix its value to ro = 2.5A.0 
With this prescription, Gersch and Rodriguez predicted a J{q, Y) that visibly overestimates 
the measured strength of the response around its maximum. This failure was later discussed 
and partially attributed to a somewhat excessively simplified approximation to the prob- 
lem.lii Nevertheless, this discrepancy seems to be eliminated by choosing a different value of 
tq. In order to show this feature, several calculations using Eqs. (|30D and (^) with different 
values of ro have been performed. In Fig. 8, results for R{q,Y) with tq equal to 2.0A, 
2.lAand 2.2Aare depicted and compared to -R(g, Y) computed with the variational p2. Even 
though the behavior of the tails of -R(g, Y) in the Gersch-Rodriguez approximation of p2 is 
different from the one of R{q., Y) with the variational two-body density matrix, the height 
and width of both peaks are in good agreement for a value of ro laying between 2.1 and 2.2A. 
Then, a proper choice of ro in the simple Gersch-Rodriguez model for p2 produces accurate 
results, provided that the height and width of the central peak are the most important 
features of the FSE broadening function. 

We have compared our results for i?(g, Y) and J(g, Y) with those obtained by SilveiH and 
Carraro and Koonin.0 Figures 9, 10 and 11 show -R(g, Y) and J(g, Y) in Gersch-Rodriguez 
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(GR), Silver (HCPT) and Carraro and Koonin (CK) theories for three values of q, 23.1 A ^, 
I5.5A-1 and 10.2A-^ The FSE function R{q,Y) is slightly different in the three theories, 
though both the height and width of the central peak are quite similar. The tails of the 
FSE broadening function show a different behavior, although they quickly vanish as | y | 
increases. Despite of the discrepancies in R{q,Y), the predicted responses are nearly the 
same at q = 23.1 A~^ and in good agreement with the experimental data. As q is lowered, 
the deficiencies of the FSE theories show up but J(g, Y) is still reasonably well described 
at g = 15.5 A~^. For the lowest q value, g = 10 A~^ (Fig. 11), the theoretical responses 
move away from experiment, and in particular do not present the small shift of the peak 
to negative Y values (see also Fig. 6). Then, even for intermediate q values, the Gersch- 
Rodriguez theory reproduces the dynamic structure function as precisely as other existing 
theories for the FSE. 

VI. SUMMARY AND CONCLUSIONS 

In this work. Final State Effects on the density response of superfiuid ^He have been 
studied in the framework of the Gersch-Rodriguez theory using a realistic description of 
the ground state of the liquid. The response is predicted as the convolution product of the 
Compton Profile JjAiY) and the FSE broadening function R{q,Y). 

Two quantities describing the ground state of the system are needed. The first one is 
the momentum distribution n{k), which completely determines Jia{Y). The second one is 
the semidiagonal two-body density matrix, which enters in the Gersch-Rodriguez form of 
Riq,Y). 

Jia{Y) has two terms, one corresponding to the non-condensate part of n{k) and another 
given by noS(Y). This splitting produces, after convoluting with R{q,Y), a total response 
which is also the sum of two terms, corresponding to the condensate and non-condensate 
contributions. The former is hnear in uq and mostly affects J(g, Y) around Y — 0. The latter 
is much less affected by FSE, although the effects are non-negligible. We have verified that 
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Gersch-Rodriguez theory gives accurate results when proper forms for the one- and two- 
body density matrices are used. A variational p2 obtained in the HNC framework accurately 
reproduces the experimental response at high g's. Furthermore, we have checked that the 
functional decomposition of p2 is very important in the calculation of R{q,Y). Simple 
models conserving the variational functional form can also produce a good estimation of the 
response. 

Our results are comparable to other calculations using more recent convolutive FSE 
theories. None of the theories correctly accounts for the observed response when q is lowered 
below about 10A~^. Further improvements could arise when higher order terms in the 
Gersch-Rodriguez cumulant expansion are considered or the time dependence of the particle 
coordinates is taken into account. 
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FIGURES 



FIG. 1. Real and imaginary parts of R{q,x) &t q = 23.1 A ^. 
FIG. 2. Functions (p{q,x) and i^{q,x) at q = 23.lA~^. 

FIG. 3. Riq,Y) atq = 23.1 A^^ (solid line) and q = 15.0 A'^ (dashed line). 

FIG. 4. Different contributions to J{q,Y) at g = 23.1 A~^. Dotted line, non-condensate term 
of JiA{y)', long-dashed line, non-condensate term of Jia{Y) after the convolution with R{q,Y); 
short-dashed line, condensate contribution once broadened by FSE; solid line, total response. 

FIG. 5. Effects of the different broadenings to the response at g = 23.1 A~^. Dotted line, 
non-condensate lA prediction; dashed line, lA broadened by FSE; solid line, total J{q, Y) including 
both FSE and IRE. 

FIG. 6. Comparison of the predicted J{q,Y) at (a) 23.1 A-^ (b) 17.9A-\ (c) IS.OA"!, and 
(d) 10.2A~^ with experimental data (points with error bars). 

FIG. 7. Detail of the central peak of the response at q = 23.1 A~^ as predicted using two 
different n(fe)'s. The solid and dashed lines correspond to njrik) and nj{k), respectively. The 
points with errorbars are the experimental data. 

FIG. 8. Comparison of i?(g, Y) calculated using the Gersch-Rodriguez p2 with different values 
of ro (solid lines) and the variational p2 (dashed line). 

FIG. 9. Comparison between Gersch-Rodriguez, Silver, and Carraro and Koonin results for 
both R{q,Y) and J{q,Y) atq = 23.1 A-^ 

FIG. 10. Comparison between Gersch-Rodriguez and Silver results for both i?(g, Y) and J(g, Y) 
at q = 15.0 A-^ 
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FIG. 11. Comparison between Gersch-Rodriguez and Carraro and Koonin results for both 
R{q,Y) and J{q,Y) at q = 10.2 A-^. 
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